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Context. In atmospheres of magnetic main-sequence stars, the diffusion of chemical elements leads to a number of observed anoma- 
lies, such as abundance spots across the stellar surface. 

Aims. The aim of this study was to derive a detailed picture of the surface abundance distribution of the magnetic chemically peculiar 
star HD 3980. 

Methods. Based on high-resolution, phase-resolved spectroscopic observations of the magnetic A-type star HD 3980 the inhomoge- 
neous surface distribution of 13 chemical elements (Li, O, Si, Ca, Cr, Mn, Fe, La, Ce, Pr, Nd, Eu, and Gd) has been reconstructed. 
The INVERS12code was used to invert the rotational variability in line profiles to elemental surface distributions. 
Results. Assuming a centered, dominantly dipolar magnetic field configuration, we find that Li, O, Mg, Pr, and Nd are mainly 
concentrated in the area of the magnetic poles and depleted in the regions around the magnetic equator. The high abundance spots of 
Si, La, Ce, Eu, and Gd are located between the magnetic poles and the magnetic equator. Except for La, which is clearly depleted in 
the area of the magnetic poles, no obvious correlation with the magnetic field has been found for these elements otherwise. Ca, Cr, 
and Fe appear enhanced along the rotational equator and the area around the magnetic poles. The intersection between the magnetic 
^ and the rotational equator constitutes an exception, especially for Ca and Cr, which are depleted in that region. 

Conclusions. No obvious correlation between the theoretically predicted abundance patterns and those determined in this study could 
1-1 be found. This can be attributed to a lack of up-to-date theoretical models, especially for rare earth elements. 

t-H Key words, stars: chemically peculiar - stars: atmospheres - stars: surface abundance structure - stars: individual: HD 3980 
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^ 1. Introduction served through rotational line profile variations of respective ele- 

L •. ments. Identifying links between abundance peculiarities, mag- 

CN Chemically peculiar (CP2) stars are main-sequence objects with netic fidd configuration; and other atmospheric parameters thus 

CD abundance peculiarities such as strong over- or underabundance helps t0 put constraints on the efficiency of diffusion processes 

(O °f certain elements compared to the atmosphere of the Sun. j n i ow - den sity plasma 
CO Compared to their "normal analogs", large magnetic fields are 

frequently detected, with characteristic field strengths of the a HD 3980 Phe or HR 183 ) is an A " l yP e star belonging 

few kG t0 S rou P °f CP2 stars as classified by Preston (1974). It 

r H Because of their rather high effective temperatures (corre- f a ^ la f fied as u an SrCrEu star b y (Bidelman & MacConnell 

X spending to spectral types B-F), relatively slow rotation, and the 1973 >- A ™ n % ^ ma ^ lc stars ° f spectral type A (Ap stars), 

H presence of strong magnetic fields, the mixing processes in at- the so " cal ed ra P ldl y oscillating (roAp) stars show fast pulsations 

^ mospheres of Ap, stars are believed to be far less efficient than in ° n ^f 1 ™f Z u Y ' ^ 

any other stars of the same spectral type. This opens a possibil- Kochukhov 2008b; Kurtz 2009). Although HD 3980 is observed 

n i , . T * * i *u u ■ i * close to the temperature range covered by roAp stars, no pulsa- 

lty for slow dynamical processes to control the chemical struc- . , t \ f ^„ * . ' \ r 

. r . , u 4-- i A-tt - i + u + - i tions have been detected so far (Weiss 1979; Martinez & Kurtz 

ture of atmospheres, where particle diffusion rules the spatial „ . , . . _ , / 

a^^u ,.r" , „ aut^L,^ . ^- . . i- r iU „ , ( , 1994; Elkm et al. 2008). A magnetic field was first observed 

distribution ol atoms ol different species, depending on the bal- . x ' ^ nnn , ^ K . ; , TTT ^™o^. • i i • 

ance between gravitational and radiative forces (Michaud 1970). in Mait f n ? l aL ( 198 °)" 3 9 80 f a vlsua 7 bl ; 

As a result, elements become inhomogeneously distributed both ^ W ". Se ?" ^ Vm§ ^ aPPafent bn S htness in V of 

9 98 mas (Weiss 1979) 

vertically and horizontally, which results in vertical abundance ' 5 v } ' 

stratification and surface abundance spots. These spots are ob- Because of its brightness and relatively fast rotation, 

HD 3980 is an ideal target for a detailed study of surface abun- 

Send offprint requests to: D. Shulyak, dance inhomogeneities. In this work we aim to reconstruct the 

e-mail: denis.shulyak@gmail.com abundance maps based on phase-resolved, high-resolution, and 
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Fig. 1. Phase coverage of the data used for Doppler imaging. The 
different panels show the phase coverage of the several wavelength re- 
gions, which changes due to the combination of spectra from various 
spectrographs and settings (see text). 



high signal-to-noise spectra obtained with different instruments. 
Dense phase and wavelength coverage allows abundance maps 
to be derived for 13 chemical elements (Li, O, Si, Ca, Cr, Mn, 
Fe, La, Ce, Pr, Nd, Eu, and Gd), which presents one of the most 
complete Doppler imaging studies to date. 



2. Observations 

Several of high-resolution spectra were obtained for studying the 
chemical structures in HD 3980 . A FEROS spectrum was ob- 
tained with the 2.2 m ESO telescope at La Silla (Chile) under 
the agreement with the Observatorio Nacional (Brazil). The Mt. 
Stromlo 74. 3 -inch telescope and the coude echelle spectrograph 
were used to obtain a series of spectra in September-October 
2001, as presented in (Polosukhina et al. 2003). Addition- 
ally, spectra acquired with UVES in 2004 (074.D-0392(A) Nes- 
vacil et al.) and 2005 (076.D-0535(A) Nesvacil et al.), HARPS 
data (074.C-0102(A) Hatzes et al., 075.C-0234(B) Hatzes et al., 
076.C-0073(A) Hatzes et al), and UVES spectra from 2006 
(077. D-0 150(A) Kurtz et al.) were included in this analysis. The 
HARPS and UVES data from 2006 were acquired via data min- 
ing and downloaded from the ESO archive 1 . 

A list of all observations is displayed in Table 1 . Since most 
of the available spectra obtained with Mt Stromlo are dedicated 
to the Li i feature at 6708 A, their spectral range is restricted to 
this wavelength region. 

Rotational phases were computed using ephemeris from 
Maitzen et al. (1980), where zero phase corresponds to the prime 
minimum in the v-filter curve of the Stromgren system (positive 
extrema of the longitudinal magnetic field): 



7X>(Prime min. in v - filter) = 



2442314.48 
0.04 + 3^9516: 



0.0003. (1) 



The relatively large rotational period error, estimated by 
Maitzen et al. (1980) from the visual examination of light curves, 
suggests significant uncertainties in the phase calculation for our 
spectroscopic data, making it problematic to relate DI maps to 
the magnetic and photometric variations. However, Elkin et 
al. (2008) shows that the phasing of the new and old magnetic 
field measurements agrees to within 0.05 of the rotational pe- 
riod. This implies a period error of « 6 x 10" 5 , which is small 
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Fig. 2. Elemental abundances of HD 3980 relative to the Sun. Vertical 
bars are attached to mapped elements and correspond to maximum and 
minimum surface abundances as shown in abundance maps (see Sect. 4 
for more details, Figs. 4, 5, and 6). Reference solar abundances are from 
Asplund et al. (2005). 



enough to be negligible for our analysis and interpretation of the 
mapping results. Our own estimates of the period accuracy led 
to a similar value of « 4 x 10" 5 . 

The resolution of the spectra ranges from 48 000 to 1 15 000, 
as shown in Table 1. For UVES spectra the resolution varies 
from 95 000 for the blue to 1 15 000 for the red arms, when us- 
ing the narrowest slit width available. Depending on the setting, 
which is different for all the UVES runs, gaps in the spectral cov- 
erage occur at different wavelengths, by combining all available 
observations the resulting phase coverage has a maximum phase 
gap of A0 = 0.2, which is visualized in Fig. 1. 

Basic data reduction was performed with pipelines of the 
corresponding instruments 2,3 . For the Mt Stromlo and FEROS 
data, the reduction was done with the Image Reduction and 
Analysis Facility (IRAF, Tody (1986)) and standard recipes. 

Additional Stromgren and Johnson photometry was used and 
is summarized in Table 3. 



3. Basic methods 

Over the course of this work two different model atmosphere 
codes were used iteratively. As a starting point, ATLAS 9 model 
atmospheres (Kurucz 1992, 1993) were computed because the 
use of opacity tables allows for very fast calculations. An initial 
abundance analysis was performed based on equivalent widths 
and individual line profile fitting in the case of blended lines. 
For further analysis LLmodels model atmospheres (Shulyak et 
al. 2004) were considered a more appropriate choice for a star 
with a rather complex abundance pattern such as for HD 3980. 
The LLmodels code was used to construct a final model atmo- 
sphere with the individual abundances presented in Fig. 2. If 
elements were found to be inhomogeneously distributed across 
the stellar surface, average homogeneous abundance values were 
implemented in the model atmosphere calculations. 

To calculate synthetic spectra the Synth3 code described by 
Kochukhov (2007) was used. The list of atomic lines used in our 
study was taken from the VALD data base (Piskunov et al. 1995; 



http ://archive . eso . org/eso/eso_archive_adp . html 



2 http ://archive .e so . org/e so/e so_archive_adp .html 

3 http://www.eso.info/sci/data-processing/software/gasgano/ 
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N. Nesvacil et al.: Multi-element Doppler imaging of the CP2 star HD 3980 
Table 1. Observations used for Doppler imaging. 



HJD 





Instrument 


Resolution 




Wavelength range [A] 


O /ICO 1 no H7C 

2 452 173.975 


n n/C/i 

0.064 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6180- 


7840 


o /ico n/i 1 cc 

2 452 174.155 


0.109 


Mt Stromlo 


nnn 

88 000 


100 


6250- 


7840 


o /ico nn nm 

2 452 170.991 


0.309 


Mt Stromlo 


nnn 

88 000 


100 


6700- 


6720 


o /ico 1 to nno 
2 452 1 /2.UU2 


n c/cc 
0.565 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


o /ico 1 to om 
2 452 1 /2.2U1 


n /c 1 c 

0.615 


Mt Stromlo 


nnn 

88 000 


1 nn 
100 


6700- 


6720 


o a co nc 1 n/i 
2 452 1 /5.194 


n oto 
0.3 11 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


o /ico 1 to 1 on 
2 452 1 /8.180 


0.128 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


o /ico 1 to non 
2 452 1 /9.080 


n c/c 
0.356 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


O /ICO 1 On A07 

2 452 180.087 


n /c 1 1 
0.611 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


/ico 1 on nco 

2 452 180.953 


0.830 


Mt Stromlo 


nnn 

88 000 


100 


6700- 


6720 


O /ICO 101 ATI 

2 452 181.0/1 


n o/cn 
0.860 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


/ico 1 00 non 

2 452 182.089 


0.117 


Mt Stromlo 


nnn 

88 000 


100 


6700- 


6720 


/ico 100 1 no 
2 452 182.192 


0.143 


Mt Stromlo 


nnn 

88 000 


100 


6700- 


6720 


a co 100 nc 1 
2 452 182.951 


n c 
0.335 


Mt Stromlo 


00 nnn 
88 000 


1 nn 
100 


6700- 


6720 


/ico ono cio 

2 453 298.518 


0.643 


T TA TTJ C n A 

UVbS 04 


11c nnn / nc nnn 

115 000/ 95 000 


350 


4800- 


6750 


/ico ion c/ci 
2 453 329.551 


n a no 
0.498 


T TA TTJ C n A 

U \hb 04 


11c nnn / nc nnn 
1 15 000 / 95 000 


ocn 
250 


4800- 


6750 


/ico oon ceo 

2 453 329.558 


0.499 


T TA TTJ C n A 

UVhS 04 


11c nnn / nc nnn 

115 000/ 95 000 


350 


4800- 


6750 


O /ICO OOO CT/I 

2 453 333.574 


0.514 


T TA TTJ O n A 

UVES 04 


11c nnn / nc nnn 

115 000 / 95 000 


300 


4800- 


6750 


/ico 000 cn 

2 453 333.577 


n c 1 c 

0.515 


T TA TTJ C n /I 

UVES 04 


11c nnn / nc nnn 

115 000/ 95 000 


onn 
300 


4800- 


6750 


O /ICO O O /I COI 

2 453 334.521 


0.754 


HARPS 


115 000 


150 


3780- 


6900 


/ico oo/i cno 

2 453 334.592 


0.772 


HARPS 


11c nnn 
115 000 


onn 
200 


3780- 


6900 


O /ICO COI T~1C\ 

2 453 581.770 


0.324 


HARPS 


115 000 


200 


3780- 


6900 


O /ICO COO 7AO 

2 453 582.798 


0.584 


T T A T> nc 


115 000 


100 


3780- 


6900 


/ico coo nn/i 
2 453 583.904 


n o/C/i 
0.864 


HARPS 


11c nnn 
1 15 000 


onn 
2UU 


3780- 


6900 


O /ICO /CO 1 Hf\H 

2 453 631.707 


0.961 


TTA^T^O nC 

UVES 05 


11c nnn / nc nnn 

115 000/ 95 000 


300 


5735 - 


9260 


O /ICO /coo /CCC 

2 453 532.555 


n om 
0.2U1 


t j\ TTj c nc 
U \hb 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


onn 
300 


5735 - 


9260 


/ico /coo /ccn 

2 453 632.659 


n ono 
0.202 


TTA^T^O nC 

UVES 05 


11c nnn / nc nnn 

115 000/ 95 000 


onn 
300 


5735 - 


9260 


O /ICO /CO/C /C/CO 

2 453 636.662 


n 1 c 
0.215 


t j\ ttj c nc 
U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


cn 
350 


5735 - 


9260 


O /ICO /CO/C /C/CC 

2 453 636.665 


n 1 c 

0.215 


tta^t^o nc 

UVES 05 


11c nnn / nc nnn 

115 000/ 95 000 


cn 

350 


5735 - 


9260 


O /ICO /CCO /CO/C 

2 453 658.626 


n TTO 

0. / /3 


t j\ ttj c nc 
U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


ocn 
250 


5735 - 


9260 


/ico /ceo /con 

2 453 658.629 


0.774 


TTA7^TZ^C nC 

UVES 05 


11c nnn / nc nnn 

115 000/ 95 000 


onn 
300 


5735 - 


9260 


O /ICO /C/CO CO/C 

2 453 662.536 


n t/co 
0. /62 


t j\ ttj c nc 
U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


onn 
300 


5735 - 


9260 


O /ICO /C/CO ceo 

2 453 665.552 


n non 
U.U2U 


tta^tztc nc 
U VE5 U5 


11c nnn / nc nnn 
1 15 UUU / 95 UUU 


onn 


5735 - 


9260 


O /ICO /C/CO CC/C 

2 455 663.556 


n non 
0.02U 


TTAft70 nC 

U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


onn 
300 


5735 - 


9260 


/ico /C/C/c T/n 

2 453 666.747 


0.828 


TTA^T^O nC 

UVES 05 


11c nnn / nc nnn 

115 000/ 95 000 


onn 
300 


5735 - 


9260 


/ico /C/ct cin 
2 453 66/. 519 


n noo 
0.023 


t j\ ttj c nc 
U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


onn 
300 


5735 - 


9260 


/ico /C/cn /cot 
2 453 669.63 / 


n ccn 
0.559 


t j\ ttj c nc 
U VE5 05 


11c nnn / nc nnn 
1 15 000 / 95 000 


onn 
300 


5735 - 


9260 


9 71 1 690 


\j. 1 oh- 


n a i\ r k j 


1 1 j yjyjyj 


900 


3780- 


6900 


2 453 712.606 


0.433 


HARPS 


115 000 


150 


3780- 


6900 


2 453 713.565 


0.676 


HARPS 


115 000 


150 


3780- 


6900 


2 453 714.622 


0.943 


HARPS 


115 000 


200 


3780- 


6900 


2 453 715.558 


0.180 


HARPS 


115 000 


200 


3780- 


6900 


2 453 716.587 


0.441 


HARPS 


115 000 


200 


3780- 


6900 


2 453 939.896 


0.952 


UVES 06 


115 000/95 000 


400 


4960- 


6990 


2 453 941.939 


0.468 


UVES 06 


115 000/95 000 


400 


4960- 


6990 


2 454 763.629 


0.407 


FEROS 08 


48 000 


300 


3560- 


9200 



Kupka et al. 1999). Further information on the final linelist used 
for reconstructing abundance maps is given in the Table 2. 

The observed line profiles were inverted with the INVERS 12 
code (Kochukhov et al. 2004). The code uses a regularized mini- 
mization algorithm and can treat many individual lines and their 
blends simultaneously. A detailed description of the software 
and some practical applications are given by, e.g., Piskunov & 
Rice (1993),Luftinger et al. (2010), and Kochukhov et al. (2004). 

4. Results 

4.1. Fundamental parameters 

The TempLogG TNG program (Kaiser 2006) was used to derive 
fundamental parameters from photometric calibrations (Moon & 
Dworetsky 1985). The available Johnson V and Stromgrem pho- 



tometry yields T Q s = 8340 K and logg = 4.3, but without any 
error estimation. 

This initial estimation was followed by an attempt to deter- 
mine r e ff and log g by fitting the Ha line with a synthetic line 
profile, which was, unfortunately, not successful owing to large 
differences between the Ha line profiles of the various spectra. 
No trend that correlates those changes to the rotational period 
was found. It is possible that these variations are caused by the 
reduction pipelines, the continuum fitting procedure, or the cur- 
vature of the echelle spectra. 

As reported in Obbrugger et al. (2008), the temperature was 
previously estimated to be 8500 K. This value was based on 
equivalent width measurements of several Fei and Fen lines, 
which also yielded log g = 4.0. Previous publications give sev- 
eral rather similar values for the temperature, e.g., 8000 ± 250 K 
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Table 2. Table of elements and spectral lines from VALD used for map- 
ping HD 3980. 



Species 




A, A 




E\aw, eV 


Blended with 


Mn ii 




-2.940 


6.128 




Mn ii 


4738.2900 


-1.876 


5.380 




Ce ii 


5518.3580 


-2.030 


0.553 




Cen 


5518.4570 


-2.000 


1.482 




Cen 


5518.4890 


-0.670 


1.155 




Cen 


5518.6030 


-2.430 


1.079 




La ii 


5880.6330 


-1.830 


0.235 


Cr 


Gdn 


5855.2150 


-1.020 


1.598 




Gd ii 


s f\r\ A CCf\f\ 

6004.5590 


-0.780 


1.659 




Fei 


6065.4820 


-1.530 


2.608 




Fei 


6230.7230 


-1.281 


2.559 


Ce 


Cm 


6176.9810 


-2.887 


4.750 




Cm 


6177.2060 


-3.093 


6.686 




Cm 


6379.7920 


-3.767 


4.497 


Fe, Gd 


Cr ii 


s r\ r\ f s c\ r\ 

6226.6380 


-3.035 


4.756 


Ce 


Eu ii 


6437.6400 


-0.320 


1.320 




Cai 


6439.0750 


0.390 


2.526 




Li i 


6707.7561 


-0.440 


0.000 


Ce, Pr, Sm 


Li i 


6707.7682 


-0.239 


0.000 




Li i 


6707.9066 


-0.965 


0.000 




Li i 


6707.9080 


-1.194 


0.000 




Li i 


6707.9187 


-0.745 


0.000 




Li i 


6707.9200 


-0.965 


0.000 




Ndra 


6145.0677 


-1.330 


0.296 


Si, Ca 


Oi 


7771.9413 


0.369 


9.146 


Cr 


Oi 


7774.1607 


0.223 


9.146 


Nd, Gd 


Oi 


7775.3904 


0.001 


9.146 




Prm 


7781.9830 


-1.210 


0.000 


Fe, Gd 


Si i 


7932.3480 


-0.468 


5.964 




Si i 


7944.0010 


-0.293 


5.984 


Cr 



(Drake et al. 2004), 8270+200 K (Kochukhov & Bagnulo 2006), 
or 8240 ± 300 K (Hubrig et al. 2007). The log g values published 
so far are 4.04 ± 0.05 from Drake et al. (2004), 4.05 ± 0.09 in 
Hubrig et al. (2007), and 4.0 ± 0.2 by Elkin et al. (2008). Conse- 
quently, T eff = 8300 ± 250 K and log g = 4.0 ± 0.2 were adopted 
for the Doppler mapping procedure in this work. 

The absolute magnitude of the star was calculated using the 
new reduction of Hipparcos data. Stellar luminosity and radius 
were computed with the bolometric correction from Landstreet 
et al. (2007). 

The value of the projected rotational velocity v sin i was de- 
termined by comparing of the synthetic spectra with observa- 
tions of magnetically insensitive iron lines. This method led to 
a v sin i= 22.5 + 2 km/s. Combination of the measured v sin i, 
period, and radius results in an inclination angle of i < 44 de- 
grees. Magnetic field measurements for this star have been re- 
ported by Maitzen et al. (1980), Hubrig et al. (2006) and Elkin 
et al. (2008). Since only measurements of the longitudinal mag- 
netic field (B\) are presently available, the magnetic field geom- 
etry is assumed to be a centered dipole (which is usually a good 
approximation in the case of magnetic A stars). Using a relation 
between the maximum and minimum of the longitudinal mag- 
netic field and inclination angle i (see Preston 1967), an angle 
J3 = 88° was therefore found between the line-of-sight and the 
rotation axis of the star. From the fit to the measurements of the 



60 6b f0 75 80 85 

P (1 

Fig. 3. Chi-square surface as a function of dipolar strength (B p ) and 
magnetic obliquity (J3). The cross marks the best-fit solution, while the 
contours correspond to the 1, 2, 3 sigma limits. The uncertainty of the 
inclination angle is taken into account implicitly. For each B p -J3 pair the 
lowest chi-square was recorded, varying inclination within its +1 sigma 
range in sin i. 



mean longitudinal field Elkin et al. (2008) obtained an amplitude 
of 1904 G. The positive magnetic pole is located near the zero 
rotational phase (light minimum). 

Because of the relatively small error in the period estimated 
by Elkin et al. (2008), we do not expect large changes in the 
position of the magnetic extrema. This conclusion also follows 
from Fig. 2 of Elkin et al. (2008), which illustrates a good match 
between the two magnetic measurements from different epochs 
(1980 and 2006) folded with the same period. Last but not 
least, the most recent magnetic measurements cited in Elkin et al. 
(2008) and our spectroscopic measurements overlap in time, thus 
justifying the use of the ephemeris from Maitzen et al. (1980). 
Finally, using the relation from Preston (1967), which combines 
the amplitude of the longitudinal magnetic field variation, in- 
clination angle, magnetic angles, and limb-darkening coefficient 
(assuming u = 0.5), we obtain the polar strength of the dipolar 
field B p « 7 kG. Taking the uncertainty of the inclination angle 
into account, the 1 sigma confidence ranges of these parameters 
are B p = 6.01 - 8.25 kG and )3 = 16° - 89°. Figure 3 illus- 
trates the analysis of the B z measurements for HD 3980, show- 
ing the chi-square surface as a function of dipolar strength (B p ) 
and magnetic obliquity (J3). 

All photometric measurements and resulting fundamental 
parameters of HD 3980 are summarized in Table 3. 



4.2. Surface abundance distributions 

In the following we give a short overview of the results for each 
individual element that has been mapped. Abundance maps are 
combined in three groups according to the distribution pattern 
and its link to the assumed magnetic field geometry. These maps 
are presented in Figs. 4, 5, and 6 (see Sect. 5). Examples of 
the theoretical fits to phase-resolved profiles of some selected 
lines are presented in Fig. 7 for illustration. When comparing 
the magnitude of over- or underabundance of some elements we 
assume reference solar abundances from Asplund et al. (2005). 
In all other cases (if not explicitly stated), when referring to the 
abundances of particular surface region, enhancements or deple- 
tion is discussed relative to the mean surface abundance. 
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Fig. 4. Abundance distributions of Li, O, Mn, Pr, and Nd at the surface of HD 3980 obtained from the lines listed in Table 2. The star is 
shown at five rotational phases. Darker areas correspond to the higher elemental abundances, the corresponding scale is given to the right of each 
panel and contours with equal abundances are plotted with steps of 1 dex. The circle and the cross indicate the position of the negative and the 
positive magnetic poles, respectively. This figure represents elements of Group 1 with enhanced spots centered on the magnetic poles. Two Li 
maps correspond (from top to bottom) to the cases of simultaneous mapping of Li i and Pr m blend at 6707 A, and additional implementation of 
Prm 7781.983 A line (see text). Abundances are given in log(^V e i/Motai) 

Lithium 



This map has the highest phase coverage ever used for Doppler 
imaging of an A-type star. Forty- six phases were available to 
reconstruct the abundance variations of Li (see Fig. 1). Prm at 
6706.703 A was varied simultaneously, because the overlap in 
wavelength occurred when Li and Pr spots were visible near the 
stellar limb. The enhancements of Li at the magnetic poles have 
an overabundance of « 5 dex compared to the Sun. The solar 



abundance of Li is « 2 dex lower than the initial cosmic abun- 
dance, (Baumann et al. 2010). 

An additional parallel implementation of Prm at 7781.983 A 
leads to a more extended concentric Li spot around the mag- 
netic poles with the same overabundance of « 5 dex (see Fig. 4, 
second plot from top). The enhancement at the negative pole 
was slightly shifted towards a higher longitude, but resulting in 
slightly worse fit to the Prm line at 6706.703 A. This could 
be due to influences of the Smu line at 6707.473 A or Cen 
6706.051 A and 6708.099 A lines that were only accounted for 



Article number, page 5 of 10 



^=0.0 



<p=Q2 



p-0.6 



^=0.8 




La 



^=0.0 a? =0.2 tf=0.4 ©=0.6 ^=0,8 , 



Eu fSgf 



p=0.6 



^=0.8 




Gd 



Fig. 5. Same as on Fig. 4, but for Si, La, Ce, Eu, and Gd. This figure represents elements of Group 2 with spots of higher abundance between the 
magnetic poles and the magnetic equator. 



as having a constant abundance over the whole stellar surface 
(see, e.g., Drake et al. 2005). 



Oxygen 

The surface abundance pattern of oxygen was reconstructed us- 
ing the triplet at 7770 A. The line profile variations are very 
pronounced and reproduced well in the area of these three O 
lines. Blends of Nd and Gd were included with constant abun- 
dance contribution. Like Li, O is concentrated in two spots at the 
magnetic poles, whereas the enhancement at the negative pole 
is also shifted towards a higher longitude. Around the magnetic 
field equator oxygen appears rather depleted. Oxygen abundance 
maps have so far been reported for the hotter Ap star e UMa 
(Rice et al. 1997) and for HD 83368 (HR 3831), with a com- 
parable temperature to HD 3980 (Kochukhov et al. 2004). It is 
interesting to note that the distribution of oxygen abundance in 
HD 3980 is completely different from what is observed in the 
two other stars, where oxygen is strongly concentrated in the re- 
gion of the magnetic equator. 



Silicon 

In order to derive the abundance pattern of Si, two Si n lines at 
7932.348 A and 7944.001 A were combined, where the weak 
Cm 7944.527 A blend was implemented by assuming a con- 
stant abundance over the rotation period. Our results indicate 
that there are several spots of overabundance around the rota- 
tional equator in the area between the magnetic poles and the 
magnetic equator. 

Calcium 

The Ca surface distribution was obtained by combining the un- 
blended Cai line at 6439.075 A with the Ndm, Can, and Sii 
blend at 6145 A. Since all elements contribute strongly to the 
blend, their distribution was determined simultaneously in this 
calculation. The resulting surface image of Ca shows an over- 
abundance in the region of the rotational equator and in the area 
of the magnetic poles of more than 2 dex compared to the Sun 
and displays approximately solar abundance around the mag- 
netic equator. 
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Fig. 6. Same as on Fig. 4, but for Ca, Cr, and Fe. This figure represents elements of Group 3 with overabundances in the area of the magnetic 
poles along the rotational equator. 



Table 3. Fundamental and atmospheric parameters of HD 3980 used for 
Doppler imaging. 



Parameter 


Value 


Reference 


Hp 


2.870 




b-y 


0.070 ± 0.007 




ml 


0.296 ± 0.007 




cl 


0.748 ± 0.022 




V 


5.71+0.005 




n, mas 


14.91+0.35 


(2) 


d, pc 


67.07 + 1.57 


(2) 


BC, mag 


0.6 + 0.07 


(3) 


Af b oi, mag 


1.64 + 0.079 




log(L/L ) 


1.24 + 0.04 




R/Re 


2.04 + 0.26 




vsini, km/s 


22.5 + 5 




U° 


60 (14-90) 




A deg 


88.1 (76-89) 




5 p ,kG 


6.9 (6.0-8.3) 




Tq&, K 


8300 + 250 




logg, cgs 


4.0 + 0.2 





Photometric measurements were taken from 1 : 
Hauck & Mermilliod (1998), 2: van Leeuwen (2007), 
3: Landstreet et al. (2007). 



Chromium 

Chromium shows a very distinct structure of an overabundance 
around the rotational equator broken by an underabundance 
around the magnetic equator. To minimize the effects of pos- 
sible vertical abundance stratification (frequently found in roAp 
stars, Shulyak et al. (2009)), two lines of the same ion with the 
similar excitation energy were used (Cm 6226.638 A and Cm 
6379.792 A). 



Manganese 

Two lines of Mnn at 4737.944 A and 4738.290 A were used. 
Due to the spectral coverage in this region, only 12 phases are 
available, but the phases are evenly spaced over the rotational 
cycle. The surface distribution of manganese is characterized by 
two patches of enhanced abundance around the magnetic poles 
and depletion at the intersections of the magnetic and rotational 
equators. 



Iron 

Two Fei lines at 6065.482 A and 6230.723 A have very similar 
excitation energy and are practically unblended. The Ce n line 
at 6231.058 A was taken into account by assuming a constant 
abundance. Iron is distributed in high-contrast spots of over- 
abundance around the rotational equator, and seems to accumu- 
late at the magnetic poles rather than at the magnetic equator. 



Lanthanum 

The most remarkable features of the surface distribution of La 
are the depletion regions around the magnetic poles. The re- 
gions with higher abundance are located in the area between the 
magnetic poles and the magnetic equator, and slightly shifted 
towards the southern hemisphere. 



Cerium 

The map resulting from the inversion of the Ce i lines at 5518 A 
shows enhanced abundances in the region between the magnetic 
poles and the magnetic equator. The area around the magnetic 
equator itself is depleted. 
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Fig. 7. Examples of observed and predicted phase-dependent profiles of selected elements computed from the respective abundance maps shown 
in Figs. 4, 5, and 6. 



Praseodymium 



The result of the inversion of Prm 7781.983 A shows that two 
spots are observed around the magnetic poles with an abundance 
increase of more than 5 dex compared to the solar value. Also the 
area of lowest abundance around the magnetic equator is more 
than 1 dex higher than for the Sun. 



We found that the combination of the Li line at 6708 A with 
the Prm 7781.983 A line results in a very similar distribution 
that shows two patches in the area of the magnetic poles, but 
their centers are shifted by +20° in longitude. The spot size, 
as well as the abundance value, agree, which is also valid for 
the depletion surrounding the magnetic equator. As mentioned 
before, the praseodymium line contributing to the Li-Pr blend at 
6708 A is rather weak. We therefore kept only the map resulting 
from the Prm 7781.983 A line. 



Neodymium 

Neodymium was - similar to Ca - derived from the combina- 
tion of the Can, Sii, and Ndm blend at 6145 A and Cai line 
at 6439.075 A; however, the contribution of Can and Sii ap- 
peared to be weak. The map shows a depletion surrounding the 
magnetic equator. Similar to Pr, two concentric areas of higher 
abundance are centered on the magnetic poles. 

Europium 

Just like Ce, the enhanced regions of Eu are situated between the 
magnetic equator and the magnetic poles. The two main spots 
are shifted towards higher longitudes relative to the magnetic 
poles. 

Gadolinium 

Gadolinium shows enhanced spots in the area of the rotational 
equator. Around the magnetic equator, several spots of lower 
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abundance are visible. The enhanced areas are located between 
the magnetic poles and the magnetic equator. This surface abun- 
dance distribution was derived by combining the two Gd n lines 
at 5855.215 A and 6004.559 A. 

5. Discussion 

5.1. Magnetic spectrum synthesis and abundance maps 

In the present Doppler imaging analysis we used a nonmagnetic 
spectrum synthesis for computating the line profiles of respective 
elements. However, with a polar magnetic field strength of up to 
7 kG and magnetic angle = 88°, the mean surface magnetic 
field modulus ranges from 4.5 kG to 5.4 kG during the rotation 
period. Such a magnetic field is strong enough to influence the 
line profiles , so that its possible impact on the abundance maps 
has to be discussed. 

It is known that the formation of the Lii 6707.473 A line is 
subject to the Paschen-Back effect (see Kochukhov 2008a; Stift 
et al. 2008), which becomes noticeable around 3 kG. However, 
we did not account for this in the present study because a) mod- 
eling the Paschen-Back effect is a non-trivial task that requires 
special software and complex techniques; b) the true geometry 
of the global magnetic field of HD 3980 is not known, i.e. its 
deviation from the simple dipole would introduce additional un- 
certainties in the line profile fitting with magnetic spectrum syn- 
thesis, and c) the star's rather high v sin /= 22.5 km/s smears out 
the line profile shape considerably. 

To estimate the impact of neglecting the magnetic field in 
spectrum synthesis on the derived abundances, especially of 
strongly overabundant elements, we used the Synth3 and Synth- 
mag codes (Kochukhov 2007) to calculate two sets of synthetic 
spectra for the lines of Fe, Cr, Eu, and Gd which are most sensi- 
tive to the magnetic field i) without magnetic field, and ii) con- 
sidering a magnetic field modulus of 5 kG. For both cases syn- 
thetic profiles were fit to observed lines in a rotation phase close 
to maximum abundance. We found that including the Zeeman 
splitting in the presence of a magnetic field leads to an abun- 
dance decrease of 0.10 - 0.15 dex for Fe and Cr. For the Eu line, 
in addition to the Zeeman splitting, isotopic splitting and hyper- 
fine structure were also taken into account in magnetic spectrum 
synthesis. A terrestrial Eu isotopic composition was adopted. 
Due to this complex pattern of the Eun A 6437 line used for 
mapping, the abundance decreased by 0.65 dex compared to the 
simplified, nonmagnetic case. The highest abundance difference 
was derived for Gd. A large Lande factor and a complex Zeeman 
pattern meant its abundance in the magnetic case had to be de- 
creased by 1.1 dex to fit observations. Except for Eu and Gd, the 
upper abundance limits of the maps computed with INVERS12 
are not very affected by neglecting the magnetic field. Lower 
abundance limits are not influenced as much for any of the lines 
used for mapping. 

Finally, it is important to understand that more accurate 
work could be done by applying a magnetic Doppler imaging 
when both distributions of the magnetic field and abundances 
are mapped simultaneously. This, of course, will become possi- 
ble once additional polarimetric observations are available. 

5.2. Surface patterns and assumed magnetic field geometry 

As noted in Sect. 4, the single-wave variations in the longitudi- 
nal magnetic field of HD 3980 can correspond to the magnetic 
dipole with the obliquity angle J3 = 88°. The true magnetic field 
configuration can, of course, be more complex. But since no 



phase-resolved observations of additional magnetic parameters 
(e.g. field modulus |B|) are available, it is currently impossible 
to determine the magnetic field topology precisely. Nevertheless, 
assuming that the magnetic field is dominated by a dipolar-like 
configuration (which is commonly found in other A-type mag- 
netic stars), we can look for a link between the global magnetic 
field geometry and elemental abundance distributions. 

The elements displayed in Fig. 4 show more or less cir- 
cular regions of higher abundance around the magnetic poles 
(Group 1). Especially around the negative magnetic pole, a shift 
in the enhanced spot towards higher longitudes is visible. Be- 
sides Pr and Nd, lighter elements such as Li, O, and Mn - an 
iron peak element - are also distributed in the same way. The 
other apparent feature of this group is lower abundance around 
the magnetic equator. 

The next set of elements (Group 2), which are displayed in 
Fig. 5 shows high abundance regions placed mainly between the 
magnetic poles and magnetic equator. Lanthanum even avoids 
the magnetic poles and shows a concentric region of underabun- 
dance around them (in contrast to, e.g., O or Pr). The region 
surrounding the magnetic equator, on the other hand, is domi- 
nated by La-depleted areas. 

The final Group 3, displayed in Fig. 6, is not as clearly de- 
fined, so it cannot be assigned to any of the other groups. Ca, 
Cr, and Fe are overabundant along the rotational equator and the 
area of the magnetic poles, but not in circular spots as in the case 
of, e.g., Li or Nd. The region of the magnetic equator is depleted 
for all those elements. 

Such a distinct separation of the mapped elements into the 
three groups described above can only be explained once the in- 
terplay between rotation, diffusion, and magnetic field is cor- 
rectly modeled in recent and/or future diffusion models. 

5.3. Comparison with diffusion models 

Assuming that the magnetic field geometry is a centered dipole, 
simple structures like rings or caps are expected, according to 
Michaud et al. (1981). Caps around the magnetic poles are found 
for Li, O, Mn, Pr, and Nd with depletion around the magnetic 
equators. The other 8 investigated elements have a more compli- 
cated structure. None of the analyzed elements show enhance- 
ments around the magnetic equator. But, e.g., Si should be over- 
abundant in areas of close to horizontal magnetic field lines, as 
stated by Alecian & Vauclair (1981), Michaud et al. (1981), or 
Vauclair et al. (1979). The enhanced spots of Si are situated be- 
tween the magnetic equator and the magnetic poles. This might 
suggest that the diffusion process is still going on. 

To explain surface abundance patterns of Ap stars, a fur- 
ther development of the diffusion theory was proposed by Ba- 
bel (1992, 1993, 1995). It is claimed that additional modifi- 
cations of the outer boundary condition by a small and inho- 
mogeneous mass outflow could considerably modify the surface 
abundance structure. Ambipolar diffusion of hydrogen in mag- 
netized plasma can also affect surface abundance patterns com- 
pared to the predictions of the initial parameter-free diffusion 
theory where no mass loss or ambipolar diffusion were consid- 
ered. In particular, for stars with a of more than 10000 K, 
a mass loss due to radiative pressure is predicted for selected 
metals (Babel 1995). Such weak metal winds could produce 
surface inhomogeneities with varying contrast between the el- 
ements. Unfortunately, only a few of the elements investigated 
in this work were the subject of diffusion calculations, which are 
so far applicable to a limited range of T q q and log g values. 
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The Group 2 of elements (e.g. Fig. 5) especially shows 
asymmetric distributions and variations of the abundance pattern 
from one element to the other. So far those cannot be explained 
by the simplified diffusion models. As stated by Kochukhov 
et al. (2004), a more complex interaction of rotation, magnetic 
field, and mass loss might be the solution for that. 

Reconstruction of the true magnetic geometry of HD 3980 by 
means of magnetic Doppler imaging (against a simplified pure 
dipolar geometry assumed in this work) is the next necessary 
step. Such an analysis will certainly provide us with deeper in- 
sight into the connection between location of the magnetic poles, 
their strengths, and distributions of chemical elements. 



6. Conclusions 

The spectrum analysis of HD 3980 was carried out in order to 
determine atmospheric parameters for Doppler imaging. By as- 
suming a centered dipolar magnetic field, the inclination of rota- 
tional axis to the magnetic field axis was determined to be close 
to 88°. Longitudinal magnetic field variation implies a dipolar 
field with a polar strength of B p « 7 kG. 

Horizontal inhomogeneities were detected and investigated 
using the Doppler imaging code INVERS12. However, this was 
done without direct incorporation of the magnetic field into the 
spectrum synthesis calculation. The resulting maps were sum- 
marized in three groups corresponding to the association of the 
surface patterns and their relation to the magnetic field. 

Lithium, O, Mn, Pd and Nd are mainly concentrated in the 
area of the magnetic poles and depleted in the regions around the 
magnetic equator. The high abundance spots of Si, La, Ce, Eu, 
and Gd are located between the magnetic poles and the magnetic 
equator. Except for La, which is clearly depleted in the area of 
the magnetic poles, no obvious other correlation with the mag- 
netic field has been found for these elements. Ca, Cr, and Fe 
show enhancements at the magnetic poles, but also along the ro- 
tational equator. The intersection between the magnetic and the 
rotational equator constitutes an exception, especially for Ca and 
Cr, which are depleted in that region. 

No obvious correlation between theoretical predictions of 
diffusion in CP stars and the abundance patterns could be found. 
This is likely attributed to a lack of up-to-date theoretical mod- 
els. For elements like O or Ca, belts of enhanced abundance 
around the magnetic equator are predicted, but not present in all 
CP stars. Unlike some CP stars, including HD 3980 analyzed in 
this paper, O is concentrated around the magnetic poles. 

Finally, to complete the investigation of the abundance struc- 
tures in the atmospheric layers of the HD 3980, further investi- 
gations will be pointed toward the accurate reconstruction of the 
surface magnetic field configuration and expected vertical abun- 
dance gradients will be investigated and modeled. 
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